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Abstract 

This article studies some numerical approximations of the homogenized matrix for stochastic linear elliptic partial 
differential equations in divergence form: — div (j4 (7,1^) Vu^{-,ijj)) = /. We focus on the case when A is a small 
perturbation of a reference periodic tensor Aper, where encodes the size of the perturbation. In this case, 
it has been theoretically shown in [6] [7| for both models considered in this article that the exact homogenized 
matrix A* possesses an expansion in powers of rj, the coefficients of which are deterministic. In practice, one cannot 
manipulate the exact terms of such an expansion. All objects are subjected to a discretization approach for the 
variables x (FEM) and cu (Monte Carlo method). Thus we need to derive a similar expansion for the approximated 
random homogenized matrix. In contrast to the expansion of the exact homogenized matrix, the expansion of the 
approximate homogenized matrix contains intrinsically random coefficients. In particular, the second order term is 
random in nature. The purpose of this work is to derive and study this expansion in function of the parameters 
of the approximation procedure (size of the truncated computational domain used, meshsize of the finite elements 
approximation) . 

1 Introduction 

The context of this work is the homogenization of stochastic hnear ehiptic partial differential equations in divergence 
form 



-div (a, V<,(-,w)) =f onV, 

u5j(-, cj) = on dV, 



(1) 



where D denotes a bounded domain of R"*, d S N* being the ambient dimension, and / e L?{T>). The matrix is 
random, symmetric, uniformly bounded and coercive, that is: 



37 > 0, e W^, £,'^Atj{x,uj)S, > 7|^P almost everywhere and ahiiost surely, 
3AI > such that |l^r,||(ioo(Rdxo))dxd < M, ^ ' 



where 7 and M do not depend on rj. Under additional stationarity hypotheses (the sense of which will be made precise 
below), it is classical that, when e — )■ 0, the random solutions of ([1]) converge in some appropriate sense to a 
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deterministic function u* solution to the fohowing problem 

r -div {A;Vu;) ^fonV, 
\ u*^0 on dV, 

where A* is the constant homogenized matrix. To compute A*, one needs first to compute some random functions, 
the so-called correctors w"^ (where p G M'*), that are solutions to random auxiliary problems, namely the corrector 
problems posed on M.'^. As a consequence, solving these corrector problems, and hence computing A*, is in general 
computationally challenging. Appropriate truncations have to be considered, and the standard numerical strategies 
lead to extremely expensive computations. In this article, we do not deal with the very general framework of stochastic 
homogenization. We rather focus on the specific case when the matrix is a small perturbation of an underlying 
periodic matrix denoted by Apcr- The parameter 77 somehow encodes the degree of randomness of A^: the smaller 77, 
the less random Ajj is. When 77 = 0, the matrix is deterministic. There are several ways to formalize the notion of a 
small perturbation in the context of stochastic homogenization (see [Tl[31|31[Slini[Zll3)- In this article, we will consider 
the following two models (for which all the details will be provided in Section 1.1 below): 

• Model 1 [Standard discrete stationary setting for stochastic homogenization]: A^(a;,cj) is stationary, and the 
small perturbation is linear (see [7]): 

Ariix,uj) = Ap^,{x) + r]Ai{x,uj) + 0{t]^), 

where the expansion holds in some appropriate functional space. 

• Model 2 [Setting introduced in 5 J: the matrix reads = Apcr (<I>^^(x,a;)), where <1>^ is referred to as a 
stochastic diffeomorphism from M.'^ to R'^ and Ap^r is a periodic matrix. An important point is that in this 
case Ajf is not necessarily stationary. The corresponding perturbative approach has been developed in [B]. The 
stochastic diffeomorphism is supposed to be a linear perturbation of the identity, namely 

$^(a;, u) = x + 77*(a;, u) + 0{rj^), 

in some appropriate functional space. 

In both settings, it has been shown that the gradients of the correctors Vw^ and the homogenized matrix A* possess 
an expansion in powers of 77: Vw^ — Vw^ + i]\/Wp + 0{ri^) and A* = Ap^^ + rjA* + 0{rf). Approximating A* by 
Apgj. + 77 one thus makes an error of the order of rf' . The functions and 77;^ are solutions to a deterministic 
and a random problem respectively. The dominant term Ap^^ in the expansion of A* is nothing but the homogenized 
matrix associated to Ap^r- The definition of A\ involves only Vtt;^ and VE (ifp). Since E (?«p) can be shown to solve a 
deterministic problem, the computation of the first order approximation Ap^^ + rjA\ only requires the solution of two 
deterministic partial differential equations. This is far less demanding that solving a single stochastic partial differential 
equation. This is the major advantage of these perturbative approaches. 

Of course, in practice, one cannot exactly solve the corrector problems that are posed on an unbounded domain. As 
will be seen below, the classical discretization approach consists in two steps: truncation of the corrector problem on a 
bounded computational domain of size N and finite elements approximation with mesh size h. This introduces two types 
of error related to the parameters N and /i, respectively. Because of the truncation, the approximated homogenized 
matrix A*^'^'^ is random. We will see that it also possesses an expansion in the variable 77, almost surely. The main 

difference between this expansion and the exact one is that, at the discrete level, the difference A^'*^'^ ~ Ap^^'^ — rjA^''^'''^ 
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is random in nature, and so is the second order error estimate. Therefore, in order to ensure the relevance of the first 
order approximation of A*-'^''^ , one needs to bound the random error r7~^(A*''''^ — A*'^,.'^ — rjAl''^'^) of second order, 

in some appropriate probability space, namely (L°°(J7))''^'' in the sequel. In the setting of Model 2, this question 
has already been addressed from a numerical point of view in The aim of the present work is to theoretically 
derive rigorous bounds on the second order error for both models and to understand how these bounds depend on the 
parameters of the discretization procedure. To this end, we first have to make precise the appropriate functional spaces 
in which the original expansions hold. We next deduce for each model in which sense the second order error terms are 
bounded. Our presentation elaborates on the previous works [51 [5]. We only provide the formalism for self-consistency 
and we refer to [6j for more details. 



1.1 Probabilistic setting 

Throughout this article, {ft, J-, P) is a probability space and we denote by E(X) ~ /^^ X{u!)clF{uj) the expectation value 
of any random variable X G L^{fl,dV). We assume that the group acts on il. We denote by {Tk)k£Zd this 

action, and assume that it preserves the measure P, that is, for all G Z'^ and all A £ T, P{TkA) = P{A). We assume 
that the action r is ergodic, that is, if A e is such that TkA = A for any k e Z"^, then ¥{A) = or 1. In addition, 
we define the following discrete notion of stationarity (see [S]): any F G ij^^. (M*^, L^{il.)) is said to be stationary if, for 
all k e Z'', 

F{x + k,io) = F{x,Tkio), almost everywhere and almost surely. (4) 

Note that this setting is a straightforward classical variant of the more commonly used (continuous) stationary 
setting for random homogenization, for which the shift r is indexed by elements of the group M.'^ and ([4]) holds for all 
k €M.'^ instead of Z'^. In our discrete setting, the ergodic theorem [TT1[T2] can be stated as follows: 

Theorem 1.1 Let F G L°° (M.'^, L^{fl)^ be a stationary random variable in the above sense. For k = {ki, k2, . . . kd) £ 
, we set |fc|oo — sup \ki\. Then 

l<i<d 

^ ^ F{x,TkUj) — > E{F{x,-)) m L'^iR'^), almost surely. 



(2N+l)d ^ ^ ' ' N 

|fc|=e<JV 

This implies that ( denoting by Q the unit cube in M."^) 

F (J^,u?j ~^^e(^J F{x,-)dx^ mL°"(R''), almost surely. 



1.2 Stochastic homogenization results 

Standard results of stochastic homogenization (see [11 [10]) apply to problem ^ when the matrix Ajj is stationary. For 
any fixed value of 77, they provide the following result. 

Theorem 1.2 We consider Model 1. Suppose that A^^ is a symmetric matrix, stationary in the sense of ([4|), uniformly 
bounded and coercive in the sense of Then the homogenized matrix A* appearing in ^ is defined by 
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where, for any p £ M.'^ , w'^ denotes the unique solution (up to the addition of a random constant) of the corrector 
problem 

— div (A,j (p + Vwp')) = on R'^, almost surely, 
Vu)^ is stationary in the sense of ^ , 



E 



= 0. 



In the sequel we will always assume that |p| = 1. 

In the case A^(a;,w) = ^por ($^"^(a;, w)), there exists an analogous result due to Blanc, Le Bris and Lions (see [51). 
Its statement requires to make precise the notion of stochastic diffeomorphism mentioned above. The map is said 
to be a stochastic diffeomorphism if it satisfies 



<I>^(-,a;) is a diffeomorphism almost surely, 
Essinf (det(V$„(x,a;))) = > 0, 

Esssup {\V^n{x,uj)\) = M' < oo, 

V^r] is stationary in the sense of (HI . 
Under these hypotheses on the following theorem gives the homogenized problem associated to 



(7) 
(8) 

(9) 

(10) 



Theorem 1.3 We consider Model 2. Suppose that Ajj — ^pcr ° '^^^ where A^^x 
bounded and coercive, and satisfies ([7]), ( 
defined by 



is a symmetric matrix, uniformly 
([9]) and (|10p . Then the homogenized matrix A*^ appearing in ([3]) is 



E 



det (V$^) 



E 



/ det (V$^) ef ^per (ej + (V$^)~' V< 



(11) 



where, for any p G , denotes the unique solution (up to the addition of a random constant) of the corrector 
problem 

-div {del (V$,,(-,w)) (V$,,(-, w))"^ ^pcr {p+ (V$,,(-, cj))"^ Vw^(-,w)j) = on W^, almost surely, 
VWp' is stationary in the sense of ([4]), ^^2) 



E 



= 0. 



1.3 Standard numerical approximation 

In practice, problems (|6]) and ()I2p are solved numerically. The first step is to introduce a truncation. Following [8], we 
approximate (|6]) by 

—div (A,,(-,w) {p + Vu;^'^(-,a;))) = on Qn almost surely, 
w''^'^ {■ , bj) is Qn — periodic, 



(13) 



where Qn — \—N~ \,N + i]"*. We consider an analogous truncated problem for (|12p . A classical finite element 
method is then used to approximate the solutions of (IT3l) . We consider a periodic triangulation Tjf^^ of the unit cell 
Q = [~^7 j]''- Replicating it, we obtain a triangulation 



Th 



iQ) 
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of R'^. We denote by V^'''^{Qn) tfie space of functions i^h defined on W^, Qjv-periodic, wliose restriction to Qjv is in a 

?Ar, and 

ifih = 0. 



'h 

^1-Lagrange finite element space built from Tfj^ = ThC^Qn, and which satisfy 



Let {(/'fc}i<fc<jv„ be a basis of V^^^ (Qn)'- 

Vpcr (Qn) = span ({0fe}i<fc<jv„) : 
where = Nv{N) is the number of degrees of freedom considered. 

In the standard case (Model 1), we define the approximated corrector w'^'^'^ as the solution to the variational 
formulation 

Find w^'^-^{-,Lo) G V^p^XQat) such that, 

yiphi^VriQN), I A^(-,t^)(p + Vu;^''''^(-,w)) -Vf^,, = almost surely, ^^^^ 

Jqn 

and the approximated homogenized matrix by 

Vl<^,J<d, K/"'^]..H = ^/^ ejA,{.,u^)(e,+Vw:f'^{-^io)). (15) 

In the second setting (Model 2), we similarly define the approximated corrector as the solution to 
Find wf''^{-,uj) e V^r'lQw) such that, for any if,, € V^,f°''(Qw), 

/ det(V$„(-,c^))(V$,,(-,c^))-^Aper(p+(V$„(-,^))"'Vu;;'^'''^(-,w)) • V^,, = almost surely, ^'^^^ 

Jqn 

and the homogenized matrix is defined by 

[^;'''^],,H=det(/Q„V$,(.,c.))"'(^/Q^det(V$,(,c.))(e, + (V$,(.,a.^ (17) 

Note that for both ([M)) and ([T6l) the solution w^'^'^ is a random field because of the truncation procedure. It follows 
that both ()15p and (|17p are random matrices. It is only in the limit N = oo that these objects become deterministic. In 
the numerical practice, one then commonly considers that the best approximation of A* is given by E (A*''''^) which 
is in turn estimated using an empirical mean, computed with standard Monte-Carlo methods. 



2 Expansion of the random matrix A^/ (Model 1) 

In this section, we consider the weakly stochastic setting of Model 1. We will turn our attention to Model 2 in Section 
3. Our main result for Model 1 is Proposition 12. II As announced above, it makes precise the behaviour of the random 
second order error in the expansion of the approximated homogenized matrix A*''^'^ , under appropriate hypotheses. 
Then, passing to the limit /i — ^ 0, we prove that this result extends to the approximated homogenized matrix when only 
truncation is taken into account. Finally, letting N go to infinity, we recover the expansion of the exact homogenized 
matrix A* derived in [7]. The functional setting described below is simple and the arguments we use in the proof 
of Proposition 2.1 are standard. Our aim is to illustrate in a simple and relevant framework that the randomness 
of the second order error in the expansions of Viu^'''''^ and A*''^'^ does not affect the validity of the approximation 

^pcr + f]^{A^'''''^). In particular, this amounts to prove that the quantity ?7~'^(A*''''^ — ^pct ~ rjA^''^'^) is bounded 
independently of /i, N and rj in some probability space. This is precisely what ensures Proposition 2.1 below. 
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2.1 Assumptions 

We suppose that the matrix A,, admits in {L°°{Q x il))'^^'^ the expansion 

A^{x,uj) = ^pcr(a;) + riAi{x,uj) + w), (18) 

where Ap^r- £ {L°°{Q)f'^, Ai e {L°°{Q x VL)f'^ and i?^ = 0{rf) in (L°°(0 x fi))''''''. This means that 

linr - Apor||(ioc(Qxa))''x'' = 0' (19) 
lim||7?-i(A„-Ap„.)-Ai||(^^(Q^^,^.,. = 0, (20) 

and there exists a deterministic constant Cn independent of rj such that, when |?7| < 1, 

ll^')ll(L-(QxO))'ix'' - '^fi' (21) 

where we recall that 

ll^ll(L~(Qxn))''xrf = Esssup \v{x,uj)\. 

{x,ui)£Qxfl 

From we know that A,, is bounded and uniformly coercive. Using ([T^ and (H]), we see that 
and, for any ^ e M'', we have 



^pcr||(ioo(Q))dxd - lim ||^r,||(ioo(QxO))''xd < M, (22) 



S,'^ Apcr{x)£, — hnr^-^A,,(a:,a;)^ > 7 |^| almost everywhere. (23) 

The matrix Apa- is thus also uniformly bounded and coercive. 

Assuming that the expansion of A,, holds in {L°° {Q x Q))'^^'^ is relevant from the point of view of modelization, 
and it somehow simplifies the proof of Proposition 2.1. As will be seen below (Remark 2.1), up to slight modifications, 
our main result (Proposition 2.1) extends to the case when the original expansion of Ajj holds in a weaker sense, namely 
when it holds almost surely with uniformly integrable bounds. Of course, in this case, the second order error in the 
expansion of A*''^''^ is bounded in a weaker sense. Lastly, it is to be mentioned that our proof does not directly extend 
to the case when the expansion of Ajj holds in spaces of the form (L°° (Q; LP{n))f'' with 0<p< +00; which models 
the idea of possibly large but rare local perturbations (see [3] for a detailed presentation of this latter model) . 

2.2 Formal expansion 

Following the method introduced in we first postulate the formal expansion of the solution to (fH|) : 

^V.h,N ^ ^0,h.N _^ ^^^1,H.N _^ ^r^,h.N^ (24) 

where Vr^''''^ — 0{rf') in some appropriate space. We will sucessively identify w^'^'^ and Wp^'^ and prove the 
validity of the expansion. Formally inserting this expression in (|14p . we obtain that the function Wp^'*'^ is independent 
of TV (we denote it by Wp'^ in the sequel) and solves 

Find w^^'' e V^P"(Q) such that. 
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and that the function Wp'^'^ is solution to 

Find G V;P"(Qw) such that, 

V^heVr'iQN), I Ap,,Wwl'''^^{-,u)-V^h+ f Ai{-,uj){p + Ww"/) -^^^=0 almost smely. ^^^^ 

Jqn Jqn 

In addition, substituting the expansions of and w"^'^'^ into (jl5p . we formaUy obtain 

A;''^'^(^) = 4''^. + + (27) 

where the terms of order zero and one are respectively defined by 

[Alt] = ^ Ape. (e, + V^«°:'') , (28) 
i.) = ^ / ef V.V.^f ^-(.,.) + ^ / ef A,(,.) (V<'' + e,) . (29) 



At 



h.N 



V \Qn\ jQf, ' ' IQwl Jq 

In the sequel, we will make precise and rigorously justify the expansions (I24|) and (I27p . 
2.3 Main result 

Our main result in this section is the following. 

Proposition 2.1 Suppose that is a symmetric matrix that satisfies (0j and is stationary in the sense of 
Suppose, in addition, that it satisfies il8\) . U9\) . \20j) . We assume that \21]) holds, namely the second order error is 
0{ri^) in {L°^{Q x Vl))'^^'^ . Then there exists a constant C independent of rj, lo, N and h, such that, for \r]\ < 1, 

?7^2 \\Vw;-''-^{;u) - Vw"/ - r^Vwl-''^^{;Lo)\\^^,^^^^.^, < Cy/\Q7\ almost surely, (30) 
where w^'^'^ , w^'^ and w^'^'^ are solutions to \14% Ii25\) and S26\) . respectively, and such that 



A*^^'^{lo) - Alt - 77At''''^(cj) < C almost surely, (31) 



where A*''''^, A^^'j and A*''*'^ are defined by (^5]) and 

Note that, as the constant C in ([50]) and PT|) is independent of ut, the expansions of Vw^'^''^ and A*''^'^ hold in 
(L°° (17;L2(Qjv)))'' and {L°°{n)f'''^, respectively. 

Proof of Proposition 2.1: Our goal is to prove estimates ([50)1 and (PT|) . To do so, following the methodology 
of [6], we begin with justifying the expansion (j24p in our discrete framework. This is the purpose of Steps 1, 2 and 
3. First we will check that Vw^''*'^ is bounded independently of 77, w, N and h and that it converges to the gradient 
WWp'^ of a deterministic function (Step 1). We will then verify that the first order error term of Ww^'^'^ , namely 
y^i),/i,Ar (j^gfined below, is bounded independently of ri,uj,N and h and that it converges to the gradient VWp'^''^ of a 
random function (Step 2). In Step 3, we prove that the second order error of S/w^'^'^ , namely Vz^''''^ defined below, 
is bounded independently of ri,ui,N and h. In other terms we prove (|30p . using bounds derived at Steps 1 and 2. In 
Step 4, remarking that the second order error term of the homogenized matrix A*''''^ depends on Vw^''*'^, Vw^''''^ 
and Vz^''*'^, we use the bounds from Steps 1, 2 and 3 to prove PT|) . 
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step 1: Our goal is first to prove that Vw^'^^^ is bounded in {^V^ {Q m))'^ , almost surely and independently of 77, w, A'^ 
and /i, and next to show that w"^'^'^ converges to a deterministic function w^'^ independent of N . Choosing iph — w^''*'^ 
as test function in p4)) and using ([2|), we obtain that 

which implies that (recall that \p\ = 1) 

\\X7„fl,h.N( ^ IIAII(L°°(Qxtj))''x'' 1^ ,1 
W^^P (•'^)|I(L^(Q„))^ ^ ^ IQ^I • 

Using we deduce that 

||Vu;^'"'^(-,c.)||(^,(^^^^, <C^|Q^|^ (32) 

where 

7 

is a constant independent of ri,uj,N and h. Thus Vw^''*'^(-, cj) is bounded uniformly in (L^((5jv))'^ independently 
from 77 and almost surely. We then remark that 

V<'^-^(-,a;) = 5^a;;_,(a.)V0., 

where = (ctp i)i<i<Af„ denotes the coordinates of it;^'''^^(-, w) in the basis of V^'^'^{Qn)- In the finite dimensional 
space V^"(QAr), all norms are equivalent and we deduce from that is bounded independently of rj in M^". 
Thus, up to extracting a subsequence, converges to a vector in M^" almost surely. Consequently, we see that 

Af„ / Ar„ \ 

V«;^-''^^(-,c^)^-^^ ^a;;,,(w)V0, = V K]<,(w)0, = Vu;°''''^(-, w) in (L2(g^.))'' almost surely, 
''^ i=i \i=i / 

where u>o^'''^(-, w) G VP"(Qw) is defined by 

i=l 

We now return to ([Tl)) which we decompose using the expansion ([T^ of A,,: for any (^/i G V^J"'' (Qat), 



JQn JQn 

(33) 

We are going to pass to the limit 77 ^ in ([55]). By definition of VWp'^'^(-, w), 

Apcr (p + Vu;2^'^'^(-, t^)) • V^h -^^ I Ape,, {p + V«;;!^'^'^(-, t^)) • V^u 
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almost surely, and 



almost surely. Using (|2T|) . we see that, when |?7| < 1, 



where we recall that is independent of 77. We thus obtain that 

lim / i?„ {p + S/wf^'^'i-^uj)) ■ V^h = 0. 
Jqn 

Passing to the limit — in (p3| . we obtain that w) satisfies 

e FP"'(Qjv), / Ap,,Vw°p'''''{;Uj) ■V^h+ f Ap,,p ■ V^h = 0. (34) 



Note that this equation has a unique solution in V^^^{Qn)- We now show that w^'^'^ = Wp''^, where w^''' is the solution 
of ([25)1 . Indeed, for any iph G V^"{Qn), we see that 

- / Ape,(Vu;^'''+p) -V I V + fc) I . 



Let 



0h := ^ iph{- + k). 



Observe that 0h is Q-periodic, and that its restriction to Q is in the Pi-Lagrange finite elements space built from T^^*^' . 
In addition, we have 



Thus 9h e VperiQ)- Hence by definition of Wp''', we obtain, for any iph G ^paiQN), 

[ Ape, (Vu;^''' + p) ■ Vifh - / Aper (Vu;°''' + p) ■ VOn = 0. 
In addition, we have w^'''' G Vp'j;r(*9) ^paiQ^)- As a consequence, Wp''* is solution to (|34l) . We conclude that 



Step 2: WG introduce the function 



V 



We want to prove that Vw^''*'^ is bounded in {L'^{Qm))'^ almost surely and independently from ri,h,N and uj, and 
that it converges to the gradient of a random function, namely w^'''''^ solution to ([25]) . Substracting equation ([Ml) 
from ([33]) . we obtain that w) g Vp^^{Qn) is such that, for any G ^'^'"^(Qjv), 

AperV<'^'^(-,w)-V^^+ / (Ai(.,a.)+rr^i?„(-,w))(p + V<'''^(-,t.))-V^;, = 0. (35) 
Choosing iph — cj) as test function in (I55t and using (^51) . we have 

\J Qn / 
< (ll^lll(L~(Qxn))'^X'^ + ll*? ^^'((•''^)||(L=o(Qj^))rfxd) 

x(||V<'"'(-,-)||(,.(q„)). + IQa^I^). (36) 

Using (PT|) . we have, when I77I < 1, 

We then deduce from ([36]), ([32]) and (l37| that, for any I77I < 1, 

< aiQjvl^ (38) 

where 

is a constant independent of r/, /i, N and w. Thus Vt'^''''^(-, w) is also bounded independently of r/ in (L^((5jv)) 
almost surely. 

As above, we deduce that there almost surely exists Wp''''^(-,w) e VpcAQN) such that, up to extracting a subse- 
quence, 

yvf^-\-,u^)^ywl-^^^^{-,u^) in {L\Q^)f. 
Using the same arguments as in Step 1, we can pass to the limit 77 — > in ([55]) . We obtain that w^''^'^ {-ju) satisfies 

Jqn Jqn 

almost surely. We thus recover ([26l) . 



d 
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Step 3: Our purpose is now to obtain a second order approximation of w^'^'^ . To this end, we define 

nnil,h,N _ p,h _ l,h,N 
- ^2 ^ V^^^{Qn)- 

Our goal is to prove that Vz^''*'^ is bounded in (^L^{Qn))'^ almost surely and independently from ri,N,h and lu. To 
do so, we will need the previous estimates (15^ and ([55]) from Steps 1 and 2 respectively. Using ([55]). p4l) and ([25]) we 
see that, for any ^pii £ V^'^''{Qn), 



We choose c^/^ = z^''*'^(-, a;) as test function in ([5^ . Using we obtain 

< ll^l|l(L-(QxO))''x'' ll^^p''''^('''^)|l(L2(Q„))d 

In view of ([21]) we have that, for any ry < 1, 

Ih ^^')(-'^)||(L-(Q„))''x'' ^ II'' ^^')||(L-(Qxn))''x'' - 

We deduce from (gO]), ([32]) and ([381) that 



Ir^p V'^-'ll(L2(Q„))'i - o^l^jvp + — (c^ + 1) IQatI" 

< (41) 

with (see ^ and ^) 

_ Pill(L~(nxQ))<'x<i / „ \ /M 1\ Cr fM 

— 



We observe that Cz is independent of N, h, rj and w. This concludes the proof of the first assertion of the proposition, 
namely the bound pOl) . 



Step 4: We now prove the second assertion in the statement of Proposition 2.1, namely the expansion of the approx- 
imated homogenized matrix. Using (|15l) and (IT51) . we have for 1 < i, j < d 

efi?,(-,^)(v</'^(-,a;) + 



1 
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Using definitions (HH), ([Ml) and we write 



1 



^pcr 



\Qn\ jQr, 

-L / eiA,{-,u)Vv:f^^A;^) 
\Qn\ Jq^ 



\Qn\ Jq 

Applying Cauchy Schwarz inequality, and the bound ([2T|). we obtain that, for any \t]\ < 1, 

1 



V 



per J jj 



V 


^*,h,N 











< 



1 



\A 



{LHQn))" 



\Al\ 



h,N 



\Qn 



(LHQn)) 



, + \Qn\). 



Thus using a.nd (??), we deduce that 







(^) 









< MCz + ||^l||(Loo(oxQ))''x'' + C'fl (1 + Cu,) . 



Recall that C^, C^, and Cz are all independent from rj, N, h and ui. The estimate ([31]) is thus proved. This concludes 
the proof of Proposition [2lTJ 



Remark 2.1 The same proof applies under weaker hypotheses. Indeed, suppose that the expansion A^i — Apcr + rjAi + 
0(772) holds m (L°°(M^)) , almost surely. This means that (|2ip is replaced by 

||77"2i?,,(-, w)||^^^^jj^j^dxd < Cr{uj) almost surely, 

where Cr is now a random variable. We suppose in addition that the random variable Cr satisfies K {\Cr\'^)'^ < Cr. 
Under this weaker assumption, we have the following estimates (compare with ([50]) and ([5T|) ]; 



with 



.l.h,N\\1 



A. 



Apci- — vAi 



Ca 



7 



< Ca, 



M 1\ Cr(uj) { M 



1 V 7 



E I 



(|mC. + \\Ai 



dxd Cv + Cr, (1 + Cw) 



\M\(L^{nycQ)ry'^ + CR{uj)'j ( ^ + ^ ) \Qn\ 
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Remark 2.2 In Proposition 2.1, the key hypothesis is (j2ip . This is the assumption from which we obtain inequali- 
ties pSI) and (|41l) (from (j36p and (j40[) respectively). Let us focus on (1401) . Looking at the first line, the term we have 
to control reads: 



\Jqn 



Using ([2T|) . we have 



1 



Qn / \JQ 

and using the hounds on 'SIw''^'^'^ previously obtained, we can conclude that Vz'^''^-'^ IS bounded in (L°° (n; L^ {Q n))) . 
Suppose now that we work within a different framework, say the original expansion of A^i holds in (^L°° (M.'^; L^(r2)))''^'^. 



Hypothesis pip is replaced by 

Esssup (e {\Rr,ix, Ol^) ' j < Cb.- 



As a consequence, we cannot expect Vz^'^'^ to be bounded in (il; ^^(Qjv)))''^'^- However, there may be a 
means to recover boundedness, but in a different space, namely {L^{Qn x fl)^'^^'^. To do so, one needs Vw^''''^ G 
{L°° {n;L°°{QN))) , which requires stronger regularity hypotheses on A^f. 

2.4 Convergence with respect to h and 

In Proposition 12. 1[ we have obtained bounds for quantities defined at the discrete level, namely after truncation and 
finite elements discretization. We now study the limit of ([30| and ([31]) when /i — > 0, and next N +oo. 



2.4.1 Convergence as h ^ Q 

First, let us define w"^'^ , and w^'^ solutions to the problems 
Find £ HI^^{Qn) such that, 

V(p e ffpi,,(Qjv), / A,(-,tj)(p + Vw;j'^(-,w)) -V^^O almost surely, ^^^^ 

Jqn 

Find w° e H^^^{Q) such that, 

JQ 

Find Wp''^{-,uj) £ H1,^^{Qim) such that, 

V(y5 e -ffp„(QAr), / AperVWp'^(-,cj) • / Ai (•, cj) (p + Vw^) • V(y5 = almost surely, 

jQw JQn 



(44) 
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respectively, where H^^^iff) denotes the closure of C'^^{Q), the space of infinitly derivable periodic functions, with 
respect to the iJ^-norm. We also define the matrices ^Jj'^, ^pcr ^i^d by 



(45) 



VI < < K'^]^^. (c) - efA,{-,u:) (e, + Vw2f{;uj)) , 

[Kor] = Aper (e, + Vwl ) , (46) 
Ur^^l (c.) = -L/ efAp„.Vu;if (,^) + -^ / efAi(.,c.)(v<+e,). (47) 



We now prove the following result, which is a direct consequence of Proposition 12. II 

Corollary 2.1 Suppose that is a symmetric matrix that satisfies and is stationary in the sense of Suppose, 
in addition, that it satisfies U8\) . f Ii20\) . We assume that i21]) holds, namely the second order error is 0{rj^) in 
{L°°{Q X f2)) ^ . Then there exists a constant C , independent of N , rj and uj such that, for \ri\ < 1, 



7?-'||Vw^'^(-,cj)-Vu;°-r/Vw^^^(-,w)||^^,^Q^^^, <C^/\Q^\ almost surely, (48) 

where w^'^, and w^'^ are solutions to ^3\ l and ^J^, respectively, and such that 

77-2 A*'^{oj) - a;,, - vAl'^ioj) < C almost surely, (49) 

where A*-'^ , A^^^. and A^''^ are defined by (|45p . (|46)) and (|47|) respectively. 
Proof : Using ((30)) , remark that 

,7-2||V<^(.,.;)-Vu;°-r/Vz«i'^(,c.)||(^,(Q^)^, < ry-^ || Vu;^''''^(., ..) ~ V^«;''^(, ^)|| (^.^q^^^. 



< V 


-2 






"V<^( 


+ V 


-2 








+ T]' 


-1 









(lHQn))" 



(50) 



where C is independent of h, N, uj and rj. Using standard properties of finite element approximations, we have that 

• lim ||Vw!J''''^(-,w) - Vw!J'^(-,'^)|Lr2,n ^^d = almost surely, 

• lim \\\/wl,-'^''^{-,uj) - Vu;^'^(-,w)|| ,,<j = almost surely, 

Passing to the limit /i — > in (|50p . we obtain The proof of follows the same line. 
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2.4.2 Convergence as N ^ oo 

We now study the limit of as — > +00. From ^ we already know that 

lim A*''^{uj) = A* almost surely, (51) 

AT— ^+00 ' ' 

where the exact homogenized matrix A* is defined by ([S])- We now turn to A^'^ whose limit when N — J> +00 is given 
by the following lemma. 

Lemma 2.1 Suppose that the matrices Aj^, Ai and Apcr o-i"^ symmetric, then the matrix A\'^ {uj) defined by (j47p 
satisfies 

lim A*{^ [lo) — A\ almost surely, (52) 
where the deterministic matrix A\ is given by 

where Wp is the unique solution (up to an additive constant) to (j43p . 
Proof : We observe that 

-.AT 



Al- 



= ^ / ef^PcrV«;if (-,0.) + -^ / efAi(,^)(v<.+e,) 
\Qn\ Jqj^ \Qn\ jQt, ^ ^ 

= (e, + VO''(Ap,rV^z;if (•,c.)+Ai(.,c.)(v<.+e,)) (53) 

because w^. is Qjv-periodic and can therefore be chosen as a test function in (|44p . Since Apct is symmetric, we have, 
using (gSl), that 

(e. + VO^VrVu'ef (-,6^)= I (V«;if (•,w))^Aper(e, + VO 



^^1. - ^ /. (^^ + V<)- A,(,.) (V< +e,) 



■JQn -^Qn 

Thus 

and we conclude using the ergodic Theorem ll.il as Ai is stationary and Vw^ is periodic. 

■ 

Using ((5T|) and (|52l) . we can pass to the limit +00 in (|49|) and we obtain that 

r^\A*^-Al^^-i^A\\<C, (54) 

where C is independent of 77 and lo (note that yl*, ^p^j. and A\ are all deterministic matrices). We thus recover the 
expansion of the exact deterministic homogenized matrix A^ as given in [7] . 
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3 Expansion of stochastic difFeomorphism (Model 2) 



We now focus on Model 2. The goal of this section is the same as that of Section 2. We prove that the random second 
order error in the expansions of the gradient of corrector Vw^''''^ and the homogenized matrix A*'''''^ is bounded 
independently of h, N and 77 in some appropriate L^lfl) space. The functional space in which the expansion of the 
original diffeomorphism holds is simple and somehow corresponds to the one considered for Model 1 in Section 2. 

3.1 Hypotheses 

In this section, we consider the Model 2 mentioned above, where the random matrix in ([1]) writes 

Ar^{x, uj) = Aper ($,7^(a;, Uj)) , (55) 

where $^ is a stochastic diffeomorphism, that satisfies conditions ([7]), ([S]), © and (llOp . The periodic matrix Ap^r is 
supposed uniformly bounded and coercive: 

37 > 0, e W^, (^Apa^{x)(^ > 7|^p almost everywhere, 
3M > such that Ppcr||(L=e(Q))dxd < M. 

This model introduced by Blanc, Le Bris and Lions in [5] is not a particular case of the standard homogenization 
setting. Following [6] we now consider a weakly stochastic case, where the diffeomorphism is close to the identity. 

More precisely, we suppose in the sequel that the following expansion holds in (C^ (IR'', L°°(51)))'' 

^r,{x,uj) = X + r]^{x,uj) + Qr,{x,LL>) with e,, = O(ry^), (57) 
where '3/ satisfies © and (|10p . This means that 

lim II*,, -Id||c.i(Rd,Loo(n))<i = 0, (58) 
l™h"'(*')-Id)-*Li(K.,Lo.(n))<* = 0, (59) 

where Id denotes the identity mapping. In addition, 0,, = 0(1]^) means that there exists a deterministic constant Ce 
independent of 77 such that, when \ri\ < 1, 

In [5], it has been shown that, under these hypotheses, the exact homogenized matrix (jlip possesses an expansion 
A* = Ap^j. + rjA\ + 0{ri^), where A*^j. and A\ are deterministic matrices that only involve solutions to deterministic 
problems posed on a bounded domain, in contrast to A* (see (|12l) and (fTTj) ). Here, as above, we will prove that the 
approximated homogenized matrix, obtained after truncation and discretization, also possesses an expansion in powers 
of 77 and that the error at second order is bounded independently of the parameters of the discretization procedure. 

3.2 Formal expansion 

As in Section 2, we first present formally our main result. We assume the following formal expansion on the approxi- 
mated corrector, solution to ([T5|: 
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where Vr^'^'^ — 0{rf') in some appropriate space. Inserting this expansion in we obtain that the function w^^^'^ 
is in fact equal to Wp''' solution to (l?5|) . The function cj) G V^por(QAf) is the unique function such that, for any 

tfh e VI\^{Qn), we have 

+ [ (div*(-,a;)Id-V*(-,a;)^)Ap„. (p + Vu;°''''^(-,tj)) •V(p,, = almost surely. (62) 

In addition, substituting Ajj and w^'^'^ by their expansions into (1171) . we formally obtain that 

A*/^^(a.) = A^.'f + T^Al^'^^'^iu) + 0(77^), (63) 
where the term of order zero is defined by ([25)1 and (1^51) , as in the first model. The term of order one is defined by 



Al 



h.N 



M = divM'(,..) + ^ /_ divvI/(,^)(e, + V<:'')^ Ape 



^ ^ (V,i;if^^(.,c.)-VvI/(.,c.)Vw;°f)'^Ape,e,. (64) 



\Qn\ Jq 

As in Section 2, we now make precise and rigorously justify the expansions (I6ip and (j63p . 
3.3 Main result 

The main result of this section is the following proposition, which is analogous to Proposition 12.1 



Proposition 3.1 Suppose that $^ satisfies Qj, 13), (Oi o-nd ilO\) . and that it is a perturbation of the identity in the 
sense of ( [57| ), Ii58\) . i59\) and I6'0)) . We suppose in addition that the symmetric periodic matrix Ap^r satisfies (j56p . 
Then there exists a constant C independent ofrj, uj, N and h, such that, for sufficiently small values ofrj, 



77-2 \\Vw;''''^{-,u) ~\7w^/- -7j\7wl'^^^{-,iu)\\^^,^^^^^, < C^/\Q7\ almost surely (65) 
where w"^'^'^ , w^'^ and Wp '^'^ are solutions to 1116]) . 1125]) and h62\]. respectively, and such that 

^ -2 - - riA];^'^{oj) < C almost surely, (66) 

where A*''''^ , A*'^^ and Al''''^ are given by ^ and ^ respectively. 

The proof, which is detailed below, follows the same lines as that of Proposition 12.11 It makes use of the following 
lemma. 



Lemma 3.1 Suppose that $^ satisfies ( |57| ), i58]) . i59]) and i60]). Then, there exist 

G (M^i°° (!}))) '^'"^ anda^ eC"^ {R'^,L°^{n)) , 

such that 

(V$^)~"^ (-,0;) = Id - riV^i{-,uj) +T,j{-,uj) almost surely, (67) 
det(V$,,)(-, cj) = 1 + r/div'i'{-,uj) + ari{-,uj) almost surely, 
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with, when rj is sufficiently small, 



(69) 



for two deterministic constants Cr and independent ofrj. In addition, there exists 770 such that, for all \ri\ < rjo , we 
have 

G M'*, Vx e M'', i|^|2<^^(V$^(x,a;))"^-(V$^(x,a;))"^C<^ICP almost surely. (70) 

Proof of Lemma 3.1: The proofs of ([ST]) . and are straightforward. We now prove ([TO]) . First, let us denote 
by A the apphcation that associates to any symmetric matrix M £ Sd(R) its ordered eigenvalues A(Af) — (Ai, • • • , A^)- 
This application is continuous on §d(K)- As a consequence, there exists a value 60 such that 



yd < 60, \M -ld\<6^ |A(M) - (1, . . . , 1)1 < i 



We infer from ^ that the following expansion holds in {C° {R'^ , L°° {n))^ : 
As a consequence, there exists a deterministic value r/o such that 



(V$^(x,w)) ^(V$^(x,w)) ^-Id 



< So, almost surely. 



Vr/ < r]o, \/x e R' 
Collecting ((TTI) and ([72]), we deduce that 

y\v\<Vo, VxeM^ |a((V$^(x,w))-^(V$^(x,w))-^) - (I,-- - ,1) 
Hence all the eigenvalues of (V$,,) (V$,,) "'^ belong to the interval [1/2,3/2]. It implies ([70 



< — , almost surely. 



(71) 



(72) 



Proof of Proposition 3.1: As announced above, the structure of the proof is the same as that of Proposition 2.1. 
First, we rigorously justify (|6T|) . Doing so, we obtain bounds on Vw^''''^, and on the error terms Vv^''^'^ and Vz^-'''^ 
that are defined below (Steps 1,2 and 3 respectively). These bounds are independent of h, N and 77. We use them at 
Step 4 to control the second order error of the expansion of the approximated homogenized matrix A*''^'^ . 

Step 1: Our goal is first to prove that Vw^'^'^ {-,0!) is bounded in {L'^{Qn))'^, almost surely. We thus show that 
^oji,N ^ the limit of 10^'^'*'^, is in fact the unique solution to ([25]) . Choosing (ph — w^''*'^(-, w) as test function in (HI 
and using (j56l) . we obtain that 



y^det(V$,)(-,w) {V%i;Uj))-'W/-^''{;L0) 



< 



det(V$^)(-,a;)|Aperp|' 



Using Lemma 13.11 and ^ we obtain 



^det(V<i>,)(-,w) iV%{;Lo))-' Vu;;'''''^(-, c.) 



> 



r],h,N 
P 



(73) 
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From dl]), we know that |det(V$^)(-,a;)| < dl{M')''. We thus have (recalhng that \p\ = 1) 

(f dct{V%){-, u)\Ap,M^y <{dKMYy M\Qn\-^. (74) 

\JQn / 

Collectmg (|73l) and (fMl) . we obtain that 

||Vu;2-"'^(-,a;)||^^,^^^^^, <C^1Q^|^ (75) 

where 

^ _ /2 (d!(M')^)^ Af 

L'lij — \ 

is independent of 77, /i, A'' and ui. Hence, Vzz;^''''^ is bounded in [QnY, independently of 77. As in the proof of 
Proposition 12 . 1 1 we deduce that there exists a function Wp''''^(-, w) g V'^'"^(QAf) such that 

V<'^'''^(-,tj) Vw°^'''^ in {L^{Qn)Y almost surely. (76) 

We now insert in equation the expansions (l67l) and (|68)) . For any yj/i e yp'j,r((3Ar), 

det(V$,(-,^))(V$^(-,a.))-^Ap„.(p+(V$,(-,L.))-iVu>;'^'''~(-,w)) -Vtp,, 
= E^^{u:)+riE^^{i^)+rfE''^{Lu) + El{i^) + F^{u), (77) 

with 

i?o'M = / Vr(p + v<'''^(-,c.)).v^„ 

£^i'(o^) = / div*(.,c.)Aper(p + Vu;f''^(-,c.)) -V^,- / VvI/(.,c.)^Ap„.(p + Vu;^'^'^(-,c.)) -V^, 

ApcrV*(-,c.)Vu;;'''''^(-,w) • V(p,„ (78) 

£;2^(c.) = / VvI/(.,c^)^Ap„.V*(-,c^)V«;^''''^(-,w)-V^^- / divvI/(., w)V*(-, l.)^ Vr (V<^''^(-, c^) + p) • 

div*(.,w)Ap„V*(-,c.)Vu;^^"'^(-,c<.) • V(^^, (79) 

£;3^(^) = / a^(-,w)Aper(p + V«;^''''^(-,w)) • V(p„+ / r,,(.,c.)^Ape,(V«;^''''^(-,c.)+p) -V.^,, 
"'Qjv JQn 

r 



Apcrr^(., ^)v<'^'''^(-, to) ■ v^h. (80) 



and F,-i{ijj) is the remainder term. Formally, Eq is of order 77", jyE'^ is of order 77, 77^-^2 is of order 77^ (because it includes 
products of two quantities of the order 77) , E^ is of order 77^ as well (it includes products of a quantity of order rj^ with 
a quantity of order rj'^) and Fjj is of order 77"^. We now prove rigorously boundedness and convergence results for each 
of the terms: Eq, E^, E2, E^ and F,,. These estimates are very useful not only at Step 1 but also at Steps 2 and 3. 
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Convergences of E^, and E'j , as 77 — !■ 0. We deduce from ([75| that 

E;1{uj)^ f A^,,{p + Vwf''''{-,u:))-V^h^ I Ap„. (p + Vz^^J^'^'^l-, w)) • V(p^ almost surely. (81) 

Jqn "'Qn 

The same argument gives 

''^^ jqn Jqn 

- AperV^'(-,w)Vu;°''''^(-,cj) • Vi^/i, almost surely (82) 
E^iuj) f V'f{-,ujfAp,,V'i>{-,uj)Vw°''''^{-,uj)-Vvh-f div*(.,c.)V*(-,w)^Ape,(V«;^^'''^(-,w)+p) • V(^„ 

-/ div^'(-,a;)AporV*(-,w)Vu;°''''^(-,tj) • V-^h, almost surely. (83) 

Bounds on E^, E2, £'3 and i^,, for sufficiently small values of 77. First we remark that 

\E1{UJ)\ < ||^pcr||(ioo(Q„))rfx. W^VhW^LHQr,)}" 

x(|QA.|H||V<^-(,.)||(,,(^^^^.)(||div*|| + llVvf ||(cO(K.,.oo(0)))-^) 

+ Ppor||(ioo(Q„))dxd ||V^'||(po(Rd^ioo(o)))''xd ||Vcp/i||(^2(Q„))<i ||V^p''''^(-:'^)||(i2(Qj^))<i : 
|£^2'(^)l < IIVr||(i^(Q„))dxd ||V(^,,||(i.(Q^))d 

X (iQjvl^ + ||Vzi;^'''^^(-,cj)||(^,(^^^^,) l|div^'||co(Rd_i»(o)) II V*||c.o(M^i-(a)) 

+ Ppcr||(ioo(Q„))dxd l|V<y5/J|(i2(Q„))d 

+ ll^*llc"(R^L-(n)) l|dlV*|lco(R<J^L°=(i2))) ' 

almost surely. Using (|75p . we deduce that there exist two constants Cei and C_b2 independent of ry, ft,, and uj such 
that, 

|£;i^(cj)| < Ce, l|V(^,,||(i2(Q„))d IQtvI^ and |£;^(cj)| < Ce^ \\^^h\\(^L^Qj,)r \Qn\K almost surely. (84) 
Second, using (|69|) . we observe that each term of E!^{uj) is bounded. As an example, the first term of E2{uj) satisfies 



a^(-,c^)Ape,(p + Vu;2^'''^(-,L.)) • 



< (7^77^ l|V¥''i|l(L2(Q„))<i Ppor||(ioo(Q„))<ixd 

x(ll^<'''(-'-)|l(L^(g„))^ + IQ^I^) 



The other terms of E^{lo) can be bounded likewise. Using (|75p . we obtain that there exists a constant independent 
of 77, /i, and w such that, for sufficiently small values of 77, we have 

\El{u)\ < C'es l|V</?/i||(i2(Q„))d IQwl^^ almost surely. (85) 

We conclude that 

E!^{uj) ^ almost surely. (86) 
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Bounding from above Fri{co) is technically more tedious, since it contains many terms. Nevertheless, each of them can 
be written under one of these two forms 

T^{;u)W/^''{;u)-^iph or / T,(.,c.)(Vz«^''''^(.,w)+p).V^,, 

where Tri{-,u}) — > strongly in {L°° {Q n))'^'^'^ almost surely. As an exemple, let us consider a term corresponding to 
the choice T^(-,cj) — ridiv'i/{-,u})ApciT^{-,uj). Using ([69)1 and d75t . we have 

/ r)diY^{;Lj)Ap,,T^(;Uj)Vw;''''^ -Vifih < ??I|div^'(-,w)Aperr„(-,w)||(^^(Q^))...C^|Qw|5 || V(p,, || (^2 (q„)). 
Qn 

< ll^pcr||(ioo(Q„))dxd ||div^'(-,cj)||^<^(Q^) Cr?7^C^,|QAr|^ || Vy^h ||(L2(Q„))d . 

This term tends to as 77 —> 0. Carrying over the same analysis for each term of we conclude that there exists 

a deterministic constant Cp independent of h, N and rj such that, for sufficiently small values of 77, we have 

< Ci^?7^ ||V(p;J|(^2(Q„))d IQatI^ almost surely. (87) 

Of course, this implies 

— > almost surely. (88) 
17— >o 

Zero order equation. We now return to equation ([77|) and pass to the limit 0. Using ([84l) . it is obvious that 

TjE'liuj) — > 0, ifE'iiuj) — > almost surely. (89) 

))— >-0 r;— >0 

Collecting dH]), dMl), (Ell) and ([Ml), we deduce that w°''''^(-, w) is a solution to JM]). We have proved (see the first 
step of the proof of Proposition 12. ip that has a unique solution which is independent of A*" and is equal to w^''^ 
the unique solution w^''^ to p5|). 



Step 2: We introduce the function 



Our aim is to prove that Vf^''''''^ is bounded in [L'^{Qm))'^ , almost surely. Then, we show that Wp'^'^ , the limit of 
^■n,h,N ^ is in fact the unique solution to (p^ . Substracting equation from ([77]) . we obtain that v"^'^'^ S ^^cAQn) 
is such that, for any G V'^^'^IQat), 



(90) 



We are now going to use the bounds on E'^, E2, E^, obtained below. Indeed, choosing Lph ~ v^''*'^(-, w) 
in ([M)) . ([SI]) . (|55|). and ([57)) . we obtain that, for sufficiently small values of 77, 

7||Vi;;'''^^(-,^)||(^,(^^^^, < (Cb, +77Ci5, +77Ci53 +772CJ.) IQjvl^ 
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where Cei, Ce2, Ces and Cp are all independent of ij, h, N and uj. As a consequence, v^'^'^{-,uj) is bounded in 
{L'^{Qn)Y almost surely. Thus, there exists a deterministic constant Cy independent of h, N and 77, such that, for 
sufBciently small values of rj, 

||Vt;^''''^(-,w)||^^2^^^^^^ < C„|Qjv|^ almost surely. 



(91) 



(92) 



Thus, there exists a function Wp-^'^ {■,uj) e Vp^^.{Qn) such that, up to extracting a subsequence, 

Vvf''^{;Uj)-^^Vwl'''^^{;uj) in{L\QN)y almost surely. 

Using dHU, dHU), dHS]), dHZl) and (HH), we can pass to the limit 77 in We conclude that w^''*'^ is the unique 

function such that, for any ifh G V^ci{Qn) 

Jom Jqm 



+ / (div*(-,w)Id- V^'(-,L^)^) Ape,(p + Vu;°''') • V^,, = 0. 



We recover the problem (|62p , the solution of which is unique. 
Step 3: We define 



C-,^) = ^ — . 



-,ri,h,N, 



Using dMl), (El), dZZl) and dm) we obtain that z"}'^^^ is such that, for any Lph £ V,f°'(QAr), 

= / Ap„Vz;-''^^(-,..).v^. + £;,^(..) + %l + :^. 

Choosing (^3^ = Zp'^'''^(-, w) in (|93| . and using (|84|) . (|85|) . and (f87|) . we obtain that, for sufficiently small values of rj, 



(93) 



7||Vz;-'^^^(.,..) 



We thus have 



|Vzp^''^-^(-,a.)| 



<C,\Qn\-\ 

where is independent of h, N, rj, and uj. Thus, inequality (j65p is proved. 



(94) 



Step 4: Our last step aims at bounding the second order error of the expansion of A*''^'^ . We insert the expan- 
sions dnZ]) and (Ull) into i^Ji). We obtain that 



\Qn\ Jq 

( 



(e, + V<''''^(-,w))^Apere, 

\- [ (e, + V<''-^(.,a.))^Ap,,e,) [ dW^i;u;) \ 

!n\ jQr, J \Qn\ Jor, 



\Ql 
1 



\Qn\ 

1 



\Q 



N\ jQ 



div*(.,c^)(e, + V<'"'^(-,c^)) A, 
(V*(-,c^)V<;'''^(-,w))^Apere, 
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where 5,^' (w) is a sum of terms of the forms 



^''(^)^ / r,(-,w)(V<:'^^^(-,w))'^Ap„.-e, and s,H-i^ / T,(.,l.) (V<:'''^(-,c.) + e,)"' Apor ■ e,, 

with s,,(aj) G M and T,, G C*^ (M"*, i°°(ri))'^^'*. For each of these (s,;,T^), there exists a constant Cs,t independent of 
/i, iV, 77, and a; such that, for sufficiently smaU values of 77, 

ll'5';llL~(f2) \\'^v\\(co{QM,L'=°in)))''^'' - ^s,T'rf'- 
Therefore, for each term of <5^'^(a;) we can write 



\Q 



N\ JQ 



< 



IQ^I ll^')ll(CO(R<',L°=(Q«)))'' 



X ||Vw;JJ''''^(-,w)||^^2(Q„))d ll^pcr||(Loo(Q„))dxd IQjV 



As a consequence, there exists a constant independent of h, N, rj and w, such that, for sufficiently small values of ry, 

\5'^'^{uj)\<Csif almost surely. (95) 



Next, we compute 



1 



\Qn\ Jq 

( ( 1 



l^wl jqn / v4m Jqn 



1 



1 



divvI/(.,c.)(V<f'^(-,w)) Ape,e, 



JQ 



pcrcj 



Using bounds (|9T|) . (|94|) and (|95p . we obtain (j66| . This concludes the proof. 



Remark 3.1 As in Section 2, we can pass to the limit h ^ in the bounds of Proposition 3.1 and prove that (|66p 
extends to the case when only truncation is taken into account. Letting N go to infinity is more difficult. Indeed, we 

need to show that lim ^lim74*''''^(a;)^ = A*, almost surely. Adapting the arguments of 181 to Model 2 will be the 

subject of a future publication. 
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